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We investigate the structure-property relationships in ABO3 perovskites exhibiting octahedral 
rotations and cooperative octahedral breathing distortions (CBD) using group theoretical methods. 
Rotations of octahedra are ubiquitous in the perovskite family, while the appearance of breath- 
ing distortions - oxygen displacement patterns that lead to approximately uniform dilation and 
contraction of the BOq octahedra - are rarer in compositions with a single, chemically unique 
B-site. The presence of a CBD relies on electronic instabilities of the B-site cations, either or- 
bital degeneracies or valence-state fluctuations, and often appear concomitant with charge order 
metal-insulator transitions or B-site cation ordering. We enumerate the structural variants obtained 
from rotational and breathing lattice modes and formulate a general Landau functional describing 
their interaction. We use this information and combine it with statistical correlation techniques to 
evaluate the role of atomic scale distortions on the critical temperatures in representative charge 
ordering nickelate and bismuthate perovskites. Our results provide new microscopic insights into 
the underlying structure-property interactions across electronic and magnetic phase boundaries, 
suggesting plausible routes to tailor the behavior of functional oxides by design. 

PACS numbers: 61.50.Ks, 31.15.xh, 71.30. +h 



I. INTRODUCTION 

Perovskite oxides with chemical formula ABO3 and 
£>-site transition metal (TM) cations exhibit a range of 
functional electronic transitions that are intimately tied 
to the structure of the fundamental building blocks: 1 (i) 
the site-symmetry of the B0 6 octahedron defined by the 
B-0 bond lengths, and (ii) the tilting of corner-connected 
octahedra. Adjacent BOq units typically fill space in per- 
ovskites through nearly rigid rotations, which produce 
deviations of the B-O-B bond angles away from the ideal 
180° found in the cubic aristotype (Pm3m symmetry); 
the rotations are described by two three-dimensional ir- 
reducible representations (irreps), and , of the 
high-symmetry structure. 2 Combinations of these lattice 
instabilities - cooperative bond length distortions and 
octahedral rotations - interact across structural phase 
transitions through elastic stresses and symmetry allowed 
coupling invariants as described within Landau theory. 

Beyond changes to crystallography, the transition 
from high temperature (high symmetry) to low tempera- 
ture (low-symmetry) can also produce electronic metal- 
insulator (MI) transitions. Perovskite oxides with B- 
site cations in t\ g e x g (<i 4 ), ^e 1 (d 7 ) and t% g ^ g (d 9 ) elec- 
tronic configurations are particularly susceptible, because 
the low-energy electronic structure is dictated by the oc- 
tahedral crystal-field split antibonding e g orbitals — the 
atomic-like d-states that are spatially directed at the 
coordinating oxygen ligands. Mi-transitions which oc- 
cur simultaneously with lattice distortions are common 
in low- dimensional materials, e.g. Peierls systems. 3 In 
three-dimensional perovskite oxides, however, the most 
familiar electronic transitions with concomitant changes 
in the B-0 bond lengths and octahedral rotations result 
from cooperative first-order Jahn- Teller effects: Tetrag- 



onal elongations of the BOq octahedra occur to remedy 
the orbital degeneracies that localized electrons encounter 
for particular B-site cation configurations. And through 
these distortions, the crystal maintains a uniform TM 
valence among all 5-sites across the transition. 4 The Jahn- 
Teller distortions are described by irreps T^, and 
i?3~ , and their interaction with octahedral rotations are 
well-established. 5 

Unusually high- valence states or e 1 electrons in delocal- 
ized band states will also produce structural distortions, 
but in most cases will preserve the uniformity of the B-0 
bonds and the octahedral crystal field of the BOq units in 
the process. 6 The B cations will readily adopt mixed va- 
lence configurations, e.g. doped perovskites manganites 7 
(containing nominally both Mn 3+ and Mn 4+ ) and also 
stoichiometric nickelates 8 (Ni 3+ ) and ferrates 9 (Fe 4+ ). 
The TM cation will not maintain an integer valence state 
(n) uniformly on all 5-sites, but rather charge dispro- 
portionate (CDP); the simplest case being into two sites 
as 

2gn+ _^ g(n-£)+ _|_ -g(n+5) + 
site 1 site 2 

where 5 is a fraction of an electron transferred between 
5-sites. Electronic correlation and on-site Coulomb re- 
pulsion effects will prefer to order the valence deviations 
S so that the inequivalent lattice sites, 1 and 2, form a pe- 
riodic arrangement, changing the translational symmetry. 
This so-called charge ordering (CO) lowers the potential 
energy of the crystal and gaps the Fermi surface. 10 

The electronic charge ordering is distinguishable 11 by 
an associated structural change in the local BOq building 
blocks (Figure 1). It appears as a B-0 bond dispro- 
portionate or "breathing distortion," which causes the 
octahedra to either dilate (site 1) or contract (site 2) ac- 
cording to the charge S transferred between sites, 12 largely 
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FIG. 1. (Color online.) Illustration of the effect of charge 
disproportionation (CDP) on the local BOq octahedral site 
equivalence. Before the CDP transition all octahedra are equiv- 
alent (left); afterwards the B06 octahedra disproportionate 
into non-equivalent sites with the simplest two-site case shown 
here: charge transfer between B-sites causes one octahedron 
to dilate (site 1) and the other octahedron to contract (site 2), 
resembling a "breathing" mode of the BOq perovskite building 
blocks. 

because of the change in ionic radii, i.e., the effective ra- 
dius of B (n ~ 5)+ is larger than B (n+5)+ . 

The magnitude of the octahedral rotations are mod- 
ified by these changes in the B-0 bond lengths. The 
extent to which cooperation between octahedral breath- 
ing and rotation instabilities is necessary to stabilizing 
charge order and MI- transit ions, however, is not as well 
understood 9,13 as the role of Jahn- Teller distortions on 
electronic transitions. 5,14 To this end, group theoretical 
methods are particularly powerful to address the interac- 
tions among the multiple octahedron-derived instabilities. 
They provide a rigorous means to evaluate the symmetry 
allowed interactions between coupled lattice degrees of 
freedom, and thus, allow us to glean insight into the mi- 
croscopic atomic structural contribution to the electronic 
CO transition in perovskite oxides. 

In this work, we enumerate the space group and order 
parameter relationships for the octahedral breathing dis- 
tortions, which are associated with irreps M±~ and of 
the aristotype cubic phase and the 15 simple octahedral 
rotation patterns available to bulk ABO3 perovskites. 
We provide a list of symmetries that perovskite oxides 
with charge order tendencies and octahedral rotations 
could adopt — not all linear combinations of the instabili- 
ties are anticipated to be symmetry allowed. 15 Using this 
information, we then illustrate how to quantify relative 
contributions of octahedral breathing distortions and ro- 
tations across phase boundaries in prototypical nickelate 
and bismuthate perovskites. This rigorous mapping of 
the unit-cell level structural distortions into a symmetry- 
adapted basis enables us to glean disentangle the role 
that the atomic structure plays in directing the macro- 
scopic electronic metal-insulator transitions. Finally, we 
show that a synergistic combination of group-theoretical 
analysis with statistical analysis makes it possible to un- 
derstand the complex interplay pervasive in mixed-metal 
perovskite oxides. 




FIG. 2. (Color online.) Schematic octahedral representations 
of the two cooperative structural breathing distortions that 
coexists with charge ordering in ABO 3 perovskites. In (a) 
the breathing distortion is two-dimensional, irrep , and 
produces a columnar ordering of the two B cations reducing the 
symmetry to PA/mmm. In (b) a three-dimensional ordering 
produces a checkerboard arrangement of the B-cations and is 
described by the irrep (Fm3m symmetry). 



II. COOPERATIVE BREATHING 
DISTORTIONS 

The octahedral breathing distortions, which create 
two unique crystallographic 5-sites, must tile in three- 
dimensions to maintain the corner-connectivity of the 
BOq framework — the defining feature of the perovskite 
crystal structure. Here we consider B-0 distortions with 
wavevectors commensurate with the Pm3m lattice peri- 
odicity, i.e., those modes which occur on the edges and 
corner of the simple cubic Brillouin zone. The two main 
cooperative breathing distortions (CBD) are illustrated 
in Figure 2. Note that the irreps, which we describe 
next, are defined using the Pm3m setting of the ABO 3 
perovskites with the B cation located at the origin. 

The first type of CBD consists of two BOq building 
blocks which are tiled to form a columnar arrangement of 
dilated and contracted octahedra [Fig. 2(a)], that splits 
the B-0 bond lengths into a doublet and quartet. 16 This 
distortion is associated with the active three-dimensional 
irrep with order parameter (a, 0, 0) and manifests as 
a zone edge k = (\, t>,0) lattice instability of the cubic 
phase. The order parameter (OP) describes a vector in 
the irrep space and corresponds to specific directions along 
which the physical distortion may be induced. For the 
irrep, the OP can have three general components (a,b,c), 
where the values a, b and c correspond to amplitudes 
of the two-dimensional breathing distortion along each 
Cartesian direction x, y and z, respectively. Here we 
consider only the case where b = c = 0, i.e. a restricted 
one-dimensional space. The possible directions for the OP 
correspond to (a, 0,0), (0, a, 0) and (0, 0, a) or columnar 
arrangements of the CBD along the x-, y- and z-directions, 
respectively, so that the octahedra distort in the same 
sense along the given direction. As a result, the symmetry 
is reduced to tetragonal, P4/mmm (space group no. 123), 
and the B-site Wyckoff position of the cubic aristotype is 
split as la — » la + lc, doubling the number of perovskite 
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TABLE I. Crystallographic data including the occupied Wyck- 
off positions (Wyck. Site) for the cooperative breathing distor- 
tions (CBD) available to ABO3 perovskites (Figure 2) in the 
absence of BOq rotations. Atom positions are given relative to 
the ideal cubic symmetry such that the relevant CBD imposed 
on the oxygen positions is indicated by A. The value of A con- 
trols the amplitude of the B—O contraction (elongation) and 
typically scales with the amount of inter-site charge transfer 
(S). The change in cell size is given relative to the pseudo-cubic 
(pc) lattice constant a pc . 
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formula units (f.u.) in the primitive cell (Table I). 

The second CBD consists of BOq octahedra which 
are tiled in a three-dimensional checkerboard arrange- 
ment [Fig. 2(b)]. The distortion is described by a one- 
dimensional irrep R± and occurs as a k = (|, |, ^) zone- 
corner lattice instability. Consequently, in the absence 
of other distortions, the CBD lifts the B-site equivalence 
while maintaining the Oh symmetry of the octahedra 
through uniform contraction and elongation of the B-0 
bonds about the trigonal axis, splitting the B-site Wyck- 
off position and forming a larger unit cell (Table I). We 
note that this type of B-site ordering is common in double 
perovskites with multiple B cations that show large ionic 
size and/or considerable oxidation state differences, e.g., 
Ba2MgW06- 17 But here, we consider only single (chemical 
species) B cations with different nominal valences. 

III. OCTAHEDRAL ROTATION AND CBD 
SPACE GROUPS 

A. Methodology 

To enumerate the allowed combinations of CBD and 
octahedral rotations in perovskite oxides we use the group 
theoretical program isotropy. 18 We follow the approach 



of Stokes et al. 19 and consider the changes in lattice 
symmetry due to the interaction of each CBD pattern 
(Figure 2) with the the 15 octahedral rotation systems 
derived from the and Rf three-dimensional irreps 
describing rotations. 2 The advantage of this approach is 
that group-sub-group relationships can be established 
between structural variants, enabling the understanding 
of the structural and electronic CO transitions within a 
first- or second-order theory. It should be noted that the 
analysis of charge order at a particular site in the crystal is 
effectively the same as analyzing the effects of cation order 
due to a order-disorder transition, because both valence 
and chemical species split a Wyckoff position identically. 
In the latter case, the irrep describes a "composition" 
mode and reflects the site occupancy. 

In this work, we enumerate those mode combinations 
with well-defined cooperative breathing distortions that 
are most likely to be observed experimentally. Although 
the rock salt B cation order has been studied previously 
in ABO 3 perovskites, the layered ordering of B cations, 
i.e. irrep M±~ , was examined only for the compositionally 
complex quadruple perovskite A4BB3O12 oxides. 20 Here, 
we report the results of the irrep on simple ABO3 
perovskites in the context of the CBD behavior, where 
cell-doubling occurs. Specifically for the irrep, we 
retain only those structures with an OP-direction (a, 0, 0) 
such that the second and third components are zero, 
unless an octahedral rotation permits it to be non-zero 
by symmetry. Furthermore, we follow the convention 
introduced in Ref. 2 and 21 and only keep structures with 
coherent rotations — those with a fixed "sense" about each 
axis. We remove structure variants from our analysis 
which would allow for a modulation in the amplitude and 
sense of the rotations about a single axis. 

B. Space Groups 

We enumerate the space groups (Table II) and associ- 
ated octahedral tilt pattern, irreps, lattice vectors, and 
origin allowed by group theory. We follow Glazer's es- 
tablished notation 22 to denote the magnitude and phase 
of octahedral tilt patterns in perovskites: The descrip- 
tion of octahedral rotations are encoded using the syntax 
a #5# c # 5 w here letters a, 6, and c indicate rotations of 
BOq units of equal or unequal magnitude about Cartesian 
y-, and z-axes. Note that in the case of equal magni- 
tude rotations about different axes, the equivalent letter 
is duplicated, e.g., aft aft aft. The superscript # can take 
three values: 0, + , or — , for no rotations, in-phase rota- 
tions (neighboring octahedra along a Cartesian axis rotate 
in the same direction), or out-of-phase rotations (adjacent 
octahedra rotate in the opposite direction), respectively. 

Table II provides the possible space group symmetries 
compatible with CBD and octahedral rotations. The 
structural data is divided into three main blocks: The 
first section contains the space group symmetries in the 
absence of CBD; the second section enumerates the sym- 
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TABLE II. Possible crystallographic space groups, octahedral rotation patterns and unit cell relationships for AB0 3 perovskites 
exhibiting rotations of octahedra given by irreps M% and R± with order parameter directions given in parentheses (r]i,r]j,rjk) 
with the CBD or Rf . The lattice vectors and origin shifts are given with respect to the high symmetry 5-atom Pm3m 
structure (B cation at the origin). 
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a In the absence of octahedral rotations, the axes are unique due to the CBD. 
b The CBD makes the z-direction unique. 

c The V denotes that this is a secondary distortion that is symmetry allowed. In this case in-phase rotations about the z-direction. The 

Glazer tilt is defined for the case b' = 0. 
d The mode acts along the body diagonal and does not create a unique axis in the crystal as in the case of the the three-dimensional 
irrep. Only the octahedral rotations lifts the axes equivalence. 



metries that result from the planar CBD with octahedral 
rotations; and the third section includes those obtained 
from the three-dimensional CBD combined with rotations. 
Structures appearing in bold in Table II correspond to 



CBD without any octahedral rotations and are given first 
at the top of block two and three. Note that for per- 
ovskites with rotations and the CBD, there is more 
than one structure possible because of the relative orien- 
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TABLE III. The change in space group which occurs when 
the different octahedral CBD are combined with the 15 simple 
octahedral tilt systems. 
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tation of the tilt pattern with respect to the columnar 
arrangement of dilated and contracted octahedra. 

Table II directly reveals the effect of the octahedral rota- 
tions and the CBD on the crystal symmetry lowering. Con- 
sider the middle block (order parameter directions M-j^, 
M3", and i?4~). The seemingly similar tilt patterns a°a°c + 
and a + a°c° - depending on the crystallographic axes they 
act upon - yield two completely different space groups, 
P2/m and Imma respectively. Furthermore, there are a 
number of space groups appearing in each of the blocks of 
Table II: C2/m appears both in the rotation only block 
(sans any CBD) and in the third block with the Rf CBD. 
Formally the crystallographic symmetries are identical 
in each case; however, physically the rotation patterns 
adopted by the crystals are different. When only 
rotations with order parameter (a, 0, b) are present in the 
crystal structure, space group C2/m allows the out-of- 
phase rotations about the y- and z-axes to be of different 
magnitude. In the presence of the three-dimensional CBD 
(Rf) only the a°b~b~ and not the a°b~c~ rotation is com- 
patible with the C2/m symmetry, whereas in the case 
of identical rotations without CBD a higher symmetry 
Imma phase results. Such restrictions imposed on the 
rotation axes equivalence are described next, and for the 
reason previously described explains why "none" appears 
as an entry in Table III. 

C. Compatible CBD and Rotation Symmetries 

In Table III, we collectively analyze how the and 
R± irreps change the space group symmetry when com- 




M+ n R+ 

a°a°a , a°a°c + , a°a° C -, a°b~b-, a~b~c-, 
a% + c~, a+a+c", a+6+c+, a%+b+ 

FIG. 3. (Color online.) Representation of the rotation pat- 
terns that are compatible with each (or both, given by the 
intersection of the) CBD without requiring a further symmetry 
reduction. 

bined with each of the 15 simple octahedral tilt systems. 
This information is schematically shown in Figure 3. Only 
one symmetries exists, PI , where the breathing distor- 
tion is geometrically compatible with the symmetry of 
the perovskite structure with octahedral rotations alone 
(a~b~c~), meaning "no change" in unit cell or transla- 
tional symmetry is required to accommodate the multiple 
distortions. Table III also reveals that the P2\/m symme- 
try is compatible with two different octahedral rotation 
patterns and the 2D M x + CBD: a°b~c- and a + b~c~. 
While PI is compatible with both and Rf irreps, 
P2\/m only exists for the irrep. 

For entries containing "none," the symmetry of the 
combined rotation pattern and CBD are incompatible, 
indicating that a further symmetry reduction would be 
required if the energetically stable phase prefers an atomic 
structure with both flavors of those distortions. This is 
made clear by examining the a~a~a~ rotation pattern, 
which corresponds to an equal amplitude of out-of-phase 
rotations about each axis, or equivalently a single out- 
of-phase rotation about the three-fold axis. The trigonal 
symmetry is incompatible with a two-dimensional distor- 
tion that would require a loss of the three-fold axis and 
therefore need at a minimum a symmetry reduction to a 
lattice with tetragonal geometry. Thus, 'none' appears 
in the column corresponding to the row with the 
a~ a~ a~ rotation pattern. Such incompatibility with the 
CBD is alleviated if the rotation pattern already 
exhibits anisotropy, e.g., a~b~b~ , yielding space group 
P2/c (Table III). 

IV. STRUCTURAL TRANSITIONS 

Although the microscopic origin for charge dispropor- 
tionation results from an electronic instability related to 
the electronic configuration of a particular metal center, 
the cooperative ordering of the CDP leads to a macro- 
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scopic bond-disproportionation. In the displacive limit, 
the two unique and nearly uniform octahedra result, di- 
lating and contracting in proportion to the magnitude of 
charge transfer. Across the electronic phase transitions, 
the breathing distortions can couple directly to octahe- 
dral rotations. (We do not consider here indirect coupling 
through a common strain component.) 

Here we describe the effect of coupling between the com- 
mon orthorhombic rotation pattern (a+b~b~) obtained 
from irreps M% © R± , which gives the six-dimensional 
order parameter with two unique directions as 

(77 + ,0,0,0,77_,77_) . 

The addition of the three-dimensional CBD Rf to the ro- 
tation pattern gives a seven-dimensional order parameter 
space, i.e., 77+ , 0, 0, 0, t?_, //_), which we contract to 
an effective four-dimensional space to obtain the P2\/c 
space group (cf. Table II). Without assuming which tran- 
sition occurs first — whether the rotations precede or follow 
the CBD — we construct a Landau free energy expansion 
about the cubic phase (/3 = t? + = rj- = 0) as 

JXA Vi) = 4? + ^ + CiA+ + C 2 p 2 r] 2 _ + Dtflrt , 

where (j) 2 and (j) 4 describe the homogeneous quadratic and 
quart ic terms for each order parameter, and C{ are the 
coefficients coupling the CBD to the in- or out-of-phase 
rotations, and D% describes the biquadratic coupling of 
the different BOq rotation "senses" in the a + b~b~ tilt 
pattern. 

The group-subgroup relationships that occur as result 
of the coupling are depicted in Figure 4. Intriguingly, 
the rotations and the CBD couple biquadratically, which 
indicates that the interactions across the transitions could 
be either cooperative or antagonistic. It is possible that 
the particular rotations could suppress CO by eliminating 
the structural octahedral breathing distortion altogether 
through the interaction term. However, in most cases, 
the rotation amplitudes are weakly modified across the 
electronic CO transition, suggesting that the strength of 
the coupling C\ is in general small. We investigate this 
in the Section V through a statistical approach. We also 
note that while the order of the coupling is important, the 
difference in temperature scales at which the rotation and 
charge ordering occurs is also important in determining 
how the structural order parameters influence each other. 
The strongest interaction occurs when the temperatures 
are similar, while if they are far apart, the two structural 
transitions will weakly couple. 



V. STRUCTURE-FUNCTIONALITY 
RELATIONSHIPS 

In this section, we apply our group theory results to 
quantitatively explore the relationship between structure 
and physical properties of experimentally known rare- 
earth nickelate and bismuthate perovskites. There is a 
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FIG. 4. Group-sub-group relationship between Pm3m and 
P2\/c. Irreps within Pm3m responsible for transitions are 
shown. Dotted line indicates that the phase transition is not 
continuous within the confines of Landau theory. 



significant interest in developing strategies - both experi- 
mentally and theoretically - to rationally control octahe- 
dral distortions through the interplay of chemical pressure, 
epitaxial strain engineering, and ultrathin superlattice 
heterostructure formation. 1,23 " 28 Although several octahe- 
dral distortion metrics, e.g., the crystallographic tolerance 
factor or bending of the B-O-B bond angle, 29 ' 30 have 
played an important role in the understanding of the elec- 
tronic and magnetic properties of perovskite oxides, 31 they 
have had limited success in materials design. Knowledge 
of quantitative structure-property octahedral distortion 
relationships are required to accelerate materials discov- 
eries. 

Distortion-mode decomposition analysis is an alter- 
native approach 5 ' 18 ' 32 (widely practiced and followed 
in the crystallography literature) to study displacive 33 
phase transitions in perovskites. It involves describing 
a distorted (low-symmetry) structure as arising from 
a (high-symmetry) parent structure with one or more 
static symmetry-breaking structural distortions. 32 In the 
undistorted parent structure, each symmetry breaking 
distort ion- mode has zero amplitude. The low-symmetry 
phase, however, will have finite amplitudes for each ir- 
rep compatible with the symmetry breaking. Said an- 
other way, the low-symmetry phase is rigorously described 
through a series expansion of the symmetry breaking 
structural modes that "freeze" into the parent structure. 
Critically, the weights or amplitudes assigned to each 
irrep are obtained according to the contribution that 
each irrep is present and the requirement that linearity 
is maintained. 

What is of particular utility in formulating quantita- 
tive relationships connecting octahedral distortions, which 
are now described mathematically, to macroscopic prop- 
erties for materials design is that each irrep carries a 
physical representation of the displacive distortions — the 
unique atomic coordinates describing various symmetry- 
adapted structural modes. The relative importance of 
these modes on properties may then be mapped by means 
of ab initio computational methods. 34 Accessibility to 
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computational methods make the distortion-mode analy- 
sis powerful, because it is possible to independently study 
various distortions (irreps) and directly assess their role 
in structural and electronic phase transition mechanisms. 
Furthermore, the distortion-mode analysis relies solely 
on crystal structural data, which enables both bulk and 
thin film materials alike to be evaluated on equal-footing. 
Such direct comparison is not possible through aggregate 
parameters such as the tolerance factor, i.e. when the 
composition is fixed, or other metrics widely followed in 
the literature. 

In the remainder of this section, we use the distortion- 
modes to form the basis for the quantitative description 
of octahedral distortions and CBD on material proper- 
ties. Using bulk i?Ni03, where R is a rare-earth element, 
and Bai-^K^BiOa perovskite compounds as prototyp- 
ical charge-ordering materials, we decompose available 
low-symmetry structural data into symmetry-adapted 
structural distortion modes. We then evaluate and cor- 
relate the amplitudes of the distortion-modes to macro- 
scopic materials behavior to uncover trends linking oc- 
tahedral distortions to the structural and physical prop- 
erties. We use the group-theory program ISODISTORT 18 
for distortion-mode decomposition analysis and R 35 ' 36 for 
the statistical analysis. For readers interested in repro- 
ducing our work, we have deposited the raw data and the 
R-script in the supplementary materials section available 
at http://link.aps.org/supplemental/XYZ. 

A. RNi0 3 Nickelates 

Rare-earth perovskite nickelates, i?Ni03, where R =Y, 
Ho, Er, Dy, Lu, Pr, or Nd, exhibit non-trivial changes in 
structure and physical properties, including sharp first- 
order temperature-driven Mi-transitions, unusual antifer- 
romagnetic order in the ground state, and site- or bond- 
centered charge disproportionation. 29,37 ' 38 At the MIT 
temperature (Tmi), the crystal symmetry lowers from 
orthorhombic Pbnm to monoclinic P2\/c symmetry 39 
(Figure 4), where the Ni cation no longer maintains a 
unique uniform valence on all sites, and disproportionates 
as 

2Ni 3+ ^ Ni^ 3 -^+ + Ni^ 3+ ^+ . 

site 1 site 2 

Moreover, the insulating ground-state displays a complex 
antiferromagnetic order (Type-E") below a Neel tempera- 
ture (T/v). Several studies 29 ' 30 ' 40 " 42 in the past have sug- 
gested the likely existence of a complex interplay between 
octhaedral rotations, transport, and magnetic properties. 
To extract deeper insight into these interrelationships, we 

(i) identify all active distortion modes in each composition, 

(ii) determine the individual amplitudes for each modes, 
and (Hi) explore the statistical correlation between in- 
dividual distort ion- modes and the physical properties, 
specifically Tmi and T/v- 



(a) Undistorted octahedra (b) M% d-type Jahn- Teller 




(c) Mg" in-phase rotation (d) out-of-phase tilting 




(e) R% a- type Jahn- Teller (f) R± out-of-phase rotation 




(g) Rc~ bending/buckling mode (h) X$~ in-phase tilting 




FIG. 5. (Color online.) Illustration of the symmetry-adapted 
orthonormal distortion modes found in the low-symmetry 
i?Ni03 and Bai-xK^BiOa perovskites. The undistorted oc- 
tahedra are shown in panel (a) for comparison. Octahedral 
tilting distortions refer to the rotation of the octahedral units 
(we use 'tilting' and 'rotation' interchangeably). The out-of- 
phase distortions are differentiated using arrows (lines and 
dashes) indicating the direction of cooperative atomic displace- 
ments, (b) M% is the d-type Jahn- Teller mode, where two 
bonds shrink and two elongate; (c) is an in-phase rotation 
mode; (d) is an out-of-phase tilting mode; (e) R% is the 
a-type Jahn- Teller mode, where four bonds contract and two 
expand, and vice versa] (f) R% is an out-of-phase rotation 
mode; (g) R£ is an out-of-phase bending mode; (h) X$ is an 
in-phase tilting mode accompanied by A cation displacements, 
which are not shown for clarity. 



1. Structure Decomposition 

We follow the procedure outlined by Campbell et 
al 18 to decompose the P2\jc monoclinic crystal struc- 
ture data, obtained from previously published diffrac- 
tion studies, 43 46 into the orthonormal symmetry-modes. 
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Diffraction data for YNi0 3 , ErNi0 3 , LuNi0 3 , and 
HoNi0 3 were measured at 295 K; DyNi0 3 at 200 K; 
NdNi0 3 at 50 K; and finally, PrNi0 3 at 10 K. T M i and 
T/v were obtained from the review article by Catalan. 30 
The P2i/c structure decomposes into eight alge- 
braically independent distortion-modes corresponding to 
the following irreps: Rt , X+, M+, M 3 + , 

and M§ . The physical representation of each irrep and 
the consequence on the octahedral framework are schemat- 
ically illustrated in Figure 5. Based on the relative am- 
plitudes of distortion-modes we obtain from isodistort, 
the dominant modes are ranked in the following order 
(from most dominant to least dominant): 

Rt > M 3 + > X+ >Rf > Rf > M+ > M+ > iJ+ 

Not surprising, the two most common octahedral 
distortion-modes that describe the a + b~b~ rotation pat- 
tern, Rt and , emerge as the dominant modes, along 
with that of X§~ : which also contains A cation displace- 
ments with relatively high amplitude. Such anti-parallel 
displacements of the A cations are established to be cor- 
related with the amplitudes of the octahedral rotations 
in perovskites. 47 There are five other irreps, namely Rr~ 
Rt , M^, and Rt with relatively low amplitudes. 

Note that the presence of the breathing mode in 
all i?Ni03 nickelates is the active mode involved in the 
symmetry lowering orthorhombic-to-monoclinic structural 
phase transition and is indicative of the CDP behavior. 

Our distortion-mode analysis points out the relatively 
weak presence of two types of Jahn- Teller distortions to 
the Ni06 units, and Rt (Figure 5). This has im- 
portant implications, because experimental observations 
based on giant oxygen 16 -0 18 isotope effect 40 indicate 
the presence of dynamic Jahn- Teller polarons; however the 
interpretation of orbital-ordering has lacked conclusive 
support from diffraction experiments. 30 More recently, 
ultra-thin films of nickelates are being increasingly investi- 
gated for rational control of orbital polarization. 48 ' 49 Our 
observation of orbital-ordering in bulk nickelates offers 
insights into the potential of tuning orbital polarization 
in epitaxial nickelate thin films. We can use the mis- 
fit strain — applied via symmetry mismatch at the thin 
film-substrate interface — to selectively enhance the contri- 
butions of Jahn- Teller modes, and Rt , in nickelates 
that already show tendencies for orbital-ordering. 50 Al- 
ternatively, orbital-ordering could emerge in nickelates, 
e.g., LaNi03 where the Jahn- Teller modes are prohibited 
by symmetry in the bulk ground-state structure, through 
quantum confinement effects as observed in the limits of 
ultra-thin films; 48,49 although in this case, strain-induced 
Jahn- Teller distortions would likely result in larger ener- 
getic penalties. 

2. Statistical Analysis 

We now evaluate the statistical correlation between 
the amplitude of the symmetry- adapted modes and the 



-1.0 -0.5 0.0 0.5 1.0 




1 1 T 



r; r; r; r; x; m; m; m; t mi t n 

FIG. 6. (Color online.) Correlation heat map capturing the 
degree of linear relationship between unit cell level structural 
distortion modes and macroscopic transition temperatures, 
Tmi, and Tn in i?Ni03 perovskites. Dark blue (red) indicates 
a strong positive (negative) correlation and white indicates 
no statistical correlation between between the two variables. 
Units for the irreps and temperatures are given in A and 
Kelvin (K), respectively. 

macroscopic Tmi, and T/v transition temperatures. We 
construct a dataset containing seven .RNiOs nickelates, 
where each rare-earth nickelate in our set is described 
using eight distortion- modes, Tmi, and T/v, resulting in 
a 7 x 10 matrix. We scaled the data by subtracting the 
mean of each column from its corresponding columns 
(this process is also known as centering) , and then divided 
the centered columns by their standard deviations. To 
evaluate the degree of linear relationship between the 
structural distortions, Tmi, and T/v, we calculate the 
sample covariance of the centered and scaled dataset as 

Cov[X] = -^f>,)(%), 
n 1 i=i 

where n = 7 is the total number of i?Ni03 compounds in 
our dataset; xi and x j are the centered and scaled column 
vectors of our data matrix X, respectively. Results from 
the covariance analysis are summarized as a correlation 
heat map in Figure 6. 

A strong positive correlation is found between modes 
capturing distortions to the Ni06 octahedra: (a pla- 
nar Jahn- Teller mode), M% , X£ , Rt, and Rt modes. 
Even though the distortion-modes are orthonormal by 



9 



(a) 



Y % 
• *Ho 

Pr 



^0.00 0.02 0.04 0.06 

r; (A) 




^0.00 0.10 0.20 0.30 

MUA) 



Dy 
Y/Ho 

» Er 



"140 150 160 

<Ni-0-Ni>(°) 



FIG. 7. (Color online.) A strong linear relationship is seen 
between (a) R% and the Neel temperature (Tn) with a correla- 
tion coefficient of 0.957 and (b) M$ and Tn with a correlation 
coefficient of 0.916. In (c) we plot the average (Ni-O-Ni) 
bond angle and Tn to illustrate that the R% and irreps 
capture the main variation in Tn for both Tmi — Tn and 
Tmi > Tn nickelates compared to the (Ni-O-Ni) bond angle. 
(Experimental data from Ref. 30.) 



construction, but when they are collectively evaluated 
for a series of RMO3 compounds our statistical analysis 
reveals that they are statistically dependent and coupled. 
These modes largely describe bond angle distortions and 
are positively correlated with the electronic transition 
Tmi- The conventional route to describe variations in 
Tmi primarily focus on tolerance factor and average Ni- 
O-Ni bond angle, whereby bending of Ni-O-Ni angle 
further from the ideal case of 180° decreases the band- 
width, promoting the insulating state over the metallic 
state. 23 ' 30 ' 51 We indeed recover this behavior, however, 
we also identify the unique displacement patterns that ge- 
ometrically sum to give the aggregate bond angle: irreps 
M2, X§ , i?^ , and cooperatively act to bend 

the Ni-O-Ni angle. These five irreps fully describe the 
Pbnm crystal structure relative to the cubic phase found 
in the metallic nickelates at high-temperature, reinforcing 
the concept that the orthorhombic distortions are largely 
responsible for the bandwidth-controlled transport behav- 
ior in nickelates, and hence prepare the electronic system 
for the MIT. 

Intriguingly, the Rf CBD mode which is the usual sig- 
nature for CO has a relatively moderate effect on Tmi, in- 
dicating that the dominate structural route to engineer the 
electronic transition may not be through isotropic bond 
length distortions. In fact, the a- type Rf Jahn- Teller 
mode and the tilting mode do not contribute signifi- 
cantly to Tmi, and we find they are anti- correlated with 
the electronic transition temperature. Our conclusion to 
this point on the Mi-transition is that the geometry of 
the oxygen framework structure in these charge-ordering 
oxides is the important atomic scale feature to tailor. 

We now shift our attention to the Type-E' antiferro- 
magnetic ordering. Unlike the transport behavior, where 
the Ni-O-Ni bond angle concept appears to be sufficient 
to explain its variability, the origin of antiferromagnetic 
ordering is much more complex. Previous high-resolution 
photoemission measurement 52 and pressure dependent 
studies 29 of Tn have established the existence of two dis- 
tinct regimes for i?Ni03 with Tmi = Tn and Tmi > Tn, 



indicative that simple tolerance factor and Ni-O-Ni su- 
perexchange interaction arguments cannot capture the 
complete physics of the system. Lee et al. used two-band 
model and Hartree-Fock theory to classify nickelates with 
Tmi — Tn as itinerant, whereas those with Tmi > Tn 
show stronger electron correlation. 53 Nonetheless, there is 
no unique descriptor known in the literature that captures 
all the variance in the Tn ordering temperature. 

Our analysis reveals the existence of a strong linear 
relationship between Tn and two irreps, R% and : the 
three dimensional Jahn- Teller mode, where in each NiOe 
octahedron four bonds contract and two elongate, and 
the out-of-phase tilting mode, respectively. The linear 
relationship is valid for both Tmi = Tn and Tmi > Tn 
nickelates (Figure 7), indicating that and con- 
tain additional information that is not captured by ei- 
ther the (Ni-O-Ni) angle or tolerance factor. The R£ 
mode is more strongly correlated with Tn than the M§ 
mode, indicating an underlying relationship between the 
long-range antiferromagnetic order and the local and sub- 
tle Jahn- Teller bond distortions. We also note that Rf 
and irreps appear in the monoclinic P2\/c crystal 
structure and are not allowed by symmetry in the Pbnm 
orthorhombic space group. 

We conjecture that with deterministic control over the 
atomic structure, e.g., enhancing R% relative to R± or 
vice versa, it could be possible to fully decouple Tmi from 
Tn, shifting the equilibirum phase boundaries dividing the 
metallic/insulating states from the paramagnetic/ordered 
regions to create non-equilibrium phase. Although it is 
compelling to extend our interpretation beyond correla- 
tion to full causation, we emphasize that first-principles 
calculations and detailed structural characterization are 
essential to establish the physical origin underlying the 
statistical associations we have identified. 



B. Bai-ccK^BiOs Bismuthates 

BaBiOs is a CO insulator which undergoes a series of 
phase transitions with temperature. At room temperature 
it has a monoclinic crystal structure with space group 
12 /m. In the monoclinic 12 /m structure, the octahedrally 
coordinated tetravalent Bi cation charge disproportionates 
into Bi( 4+(5 ) + and Bi^ 4_5 ^ + . Previous optical conductivity 
measurements 54 attribute the appearance of an electronic 
gap in the stoichiometric 12 /m phase to structural CBD 
and octahedra rotations. 

Chemical doping of BaBiOs with potassium by random 
substitution on the Ba sites, Bai-^K^BiOa, produces a 
sequence of structural and electronic phase transitions. 
Increasing the potassium concentration leads to crystal 
structures of higher symmetry: An orthorhombic Ibmm 
phase is stable over the range 0.12 < x < 0.37. At x « 
0.37, a tetragonal 14/ mem phase is stable and for 0.37 < 
x < 0.53 BaBi0 3 transforms to Pm3m cubic. 55 ' 56 The 
electronic transport properties also evolve concomitant 
with these structural changes: The system transforms 
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FIG. 8. (Color online.) Relationship between crystal symmetry and structural distortions in Bai-xK^BiOs as a function of 
A site doping. The number of symmetry unique structural modes decreases with increasing K concentration. Pure BaBiOs 
is described by four modes, R± , R% , R^, and K£ . As the concentration of K increases, the Jahn- Teller R% mode disappears 
rapidly, followed by loss of the three-dimensional CBD Rf. The bar graph (inset) shows quantitatively in units of angstroms 
the change in amplitude of the structural modes observed in each space group. 



from a robust insulating state in the monoclinic phase to 
semiconducting in the orthorhombic phase, and finally to 
a superconducting state in the tetragonal and undistorted 
cubic phases. 

Our focus here is to systematically track the struc- 
tural phase transition as a function of K doping, with 
emphasis placed on the role of the symmetry-adapted 
modes and CBD-rotation symmetries. We decompose 
the low-symmetry crystal structures of Bai_^K x Bi03 
across the phase boundaries using five chemical com- 
positions: BaBiOs, Bao.gKo.iBiOs, Bao.8Ko.2Bi03, 
Ba .63Ko.37Bi0 3 , and Ba .6oKo.4oBi0 3 , for which high- 
resolution diffraction data exists. 55,57 

The results from our analysis are summarized schemat- 
ically in Figure 8. The 12 /m structure is described by 
four irreps Rf , R% , R± , and Rf , which capture the CBD, 
a Jahn- Teller distortion, and out-of-phase octahedral ro- 
tations, and out-of-phase bond stretching, respectively. 
Although the 12 jm monoclinic symmetry is maintained 
at small doping, as with Bao.gKo.iBiOs, our mode de- 
composition analysis reveals that the amplitude of the 
R% mode vanishes completely. The relative amplitudes 
of the Rf breathing mode and out-of-phase octahedral 
rotation modes also decrease, whereas the amplitude of 
the out-of-phase bending distortion increases. 

With further increase in K concentration, for example in 
Bao.8Ko.2Bi03, the crystal structure becomes orthorhom- 
bic Ibmm. Interestingly, even in the absence of the CBD 



distortion, which we find has zero amplitude at x = 0.2, 
the electronic structure remains insulating. At the same 
time, the amplitude of the R§ mode continues to increase 
in the presence of out-of-phase octahedral rotations R± . 

Upon further increase in K, the Rr~ mode disappears 
completely, yielding tetragonal I A/ mem Bao.63Ko.37Bi03 
with only the R^ octahedral rotation mode at reduced am- 
plitude. Across this structural transition, the compound 
becomes metallic and superconducting 55 with a critical 
transition temperature of 27 K. Our analysis suggests 
that the subtle Rf distortion may play a more important 
role in the transport properties of Bai-^K^BiOs than the 
the commonly believed Rf CBD which vanishes immedi- 
ately at small potassium doping. Finally, for x~0.4, the 
superconducting phase is also stable in the cubic Pm3m 
structure without any octahedral rotations or distortions. 

Our group theoretical approach applied to the 
Bai-^K^BiOs system, reveals systematic trends in the 
irreps with doping across the structural phase bound- 
aries. One of the outstanding questions regarding the 
electronic properties of Bai-^K^BiOa is the stability of 
the wide semiconducting or insulating region with doping. 
Recently Franchini et al. performed density functional 
calculations to understand this feature of the phase dia- 
gram, attributing the insulating state to the formation 
of polaronic defects, 56 which are consistent with spectro- 
scopic measurements. 58 Here, we suggest the existence 
of a plausible relationship between the stability of the 
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polaronic defects and the persistence of the R$ distortion 
in Bao.8Ko.2Bi03, which eventually delays the onset of 
superconducting phase in Bai-^K^BiOs perovskite. 

VI. SUMMARY 

We applied group theoretical methods to evaluate the in- 
terplay between octahedral rotations and electronic charge 
ordering, which we parameterized structurally as octahe- 
dral "breathing" or cooperative bond length distortions 
through irreducible representations. We enumerated the 
possible space groups available to perovskites from the 
intersection of these distortions. This crystallographic 
data should prove useful for experimental structure refine- 
ments or serve as a well-defined set of symmetry-unique 
structures for which the stability of transition metal com- 
pounds susceptible to charge ordering may be evaluated 
with first-principles total energy methods. 

With this information, we combined structural mode 
decomposition techniques with statistical methods to ex- 
tract structure-property relationships from available ex- 
perimental diffraction data on representative charge or- 
dering oxides. The crystal structures of nickelates and bis- 
muthates were decomposed in terms of symmetry adapted 
distortion-modes and the evolution in the amplitude of 
each local structure mode tracked as a function of an 
external chemical parameter. We illustrated that this 
alternative set of descriptors provides a useful construct 
beyond the traditional tolerance factor paradigm found in 
perovskites to understand the atomic scale origin of physi- 
cal properties, specifically how unit cell level modifications 
correlate with macroscopic functionality. 

Our statistical analysis uncovered previously unappre- 
ciated relationships that may be harnessed for electronic 
structure by design, some of which are being carried out 
within our group and will be reported later. We emphasize 
that the application of statistics to the decompositions 
does not require solely experimental data; the analysis 
may be performed with computationally obtained data 
alone. Importantly, the relationships established with 



these techniques may be cross- validated by the construc- 
tion of hybrid data sets, which combine theoretical results 
with experiment data, making it possible to extract and 
validate new insight into the material physics of oxides 
with correlated electrons. We anticipate that this ap- 
proach will spawn a number of additional studies in diverse 
crystal classes since it is immediately generalizable: the 
synergy of applied group theoretical methods with statis- 
tical analysis and subsequent first-principles calculations 
provides a platform to achieve rational structure-driven 
design of complex materials. 
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Appendix A 



TABLE IV. Abbreviations appearing throughout the text. 



CBD Cooperative octahedral breathing distortion 

TM Transition metal (referring to a cation) 

MIT Metal— to— insulator transition 

CDP Charge disproportionate 

Tmi Metal— to— insulator transition temperature 

Tn Neel temperature 

DFT Density Functional Theory 

Tc Superconducting transition temperature 
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